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ABSTRACT 

We have spectroscopically confirmed a brown dwarf mass companion to the hydrogen atmo- 
sphere white dwarf NLTT 5306. The white dwarf's atmospheric parameters were measured 
using Sloan Digital Sky Survey and X-Shooter spectroscopy as T^s = 7756 ± 35K and 
log(g)= 7.68 ± 0.08, giving a mass for the primary of AfwD = 0.44 ± 0.04 Mq at a distance 
of 71 ± 4pc with a cooling age of 710±50Myr The existence of the brown dwarf secondary 
was confirmed through the near-infrared arm of the X-Shooter data and a spectral type of dL4- 
dL7 was estimated using standard spectral indices. Combined radial velocity measurements 
from the Sloan Digital Sky Survey, X-Shooter and the Hobby-Eberly Telescope's High Reso- 
lution Spectrograph of the white dwarf gives a minimum mass of 56±3 Afj„p for the secondary, 
confirming the substellar nature. The period of the binary was measured as 101 .88 ±0.02 mins 
using both the radial velocity data and i'-band variability detected with the INT. This variabil- 
ity indicates 'day' side heating of the brown dwarf companion. We also observe Ha emission 
in our higher resolution data in phase with the white dwarf radial velocity, indicating this sys- 
tem is in a low level of accretion, most likely via a stellar wind. This system represents the 
shortest period white dwarf + brown dwarf binary and the secondary has survived a stage of 
common envelope evolution, much like its longer period counterpart, WD 0137— 349. Both 
systems likely represent bona-fide progenitors of cataclysmic variables with a low mass white 
dwarf and a brown dwarf donor 
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1 INTRODUCTION 

As the descendants of high mass ratio binaries, brown dwarf 
(BD) companions to white dwarfs (WDs) enable investigation 
of one extreme of binary formation and evolution, including 
the known deficit of BD companions to main sequence stars 
jMcCarthv & Zuckermanll2004lGrether & Lineweaverll2006l) . The 
BDs can be directly detected relatively easily at all separations, 
since they dominate the spectral energy distribution at near- to mid- 
infrared (IR) wavelengths, in stark contrast to main sequence + 
substellar pairs. 

As d etached companion s to WDs, BDs are rare 
(~ 0.5%. ISteele etd] boill : iGirven et al.l l201lh . Only a 
handful of such systems have thus f ar be en spectroscopi- 
cally confirmed e.g. GDI 65 (DA +dL4. [ Bec klin & Zuckerman 
19881). GD 1 400 (DA+dL6-7, iFarihi feChristopher. ,2004t 



DobbieetalJ l2005h. WD0 137-349 (DA+dL8, iMaxted etal 



2004 iBurleigh et al.l [200^ , PHL5038 (DA+dL8 nSteele et al 



|2009|) . and LSPM 1459+0857 (DA+T4.5, lDav- Jones et al]|201lh . 
GD165, PHL5038 and LSPM 1459+0857 can be classed as 
widely orbiting with projected separations of 120AU, 55 AU 
and 16500-26500 AU respectively. WD 0137-349 and GD 1400 
have muc h shorter o rbital periods of 116 mins and ~ lOhrs 
(Burleigh e t alj I20T1I) respectively, and have both undergone a 
common envelope evolution. 

These two distinct populations are thought to be the out- 
come of stellar evolution; the wide pairs where the secondary has 
migrated outwards due to the mass los s of the WDs progenitor 
jparihi et al.ll2006l;lNordhaus et al.llioiol) , and the close systems in 
which the secondary has survived a stage of common envelope evo- 
lution and ma y eventually lead to the formation of a cataclysmic 
variable (CV) ( |Politandl2004ah . In these close binaries, the BD is 
expected to be irradiated by the WD's high UV flux, leading to 
substantial differences in the 'day' and 'night' side hemispheres. 
These systems can additionally be used for testing models of irra- 
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Table 1. SDSS and UKIDSS magnitudes of NLTT5306. 



Band 


Magnitude 


u' 


17.51 ±0.01 


9' 


17.03 ±0.00 


r' 


16.95 ±0.00 


i' 


16.96 ±0.01 


z' 


17.00 ±0.01 


Y 


16.49 ±0.01 


J 


16.24 ±0.01 


H 


15.86 ±0.01 


K 


15.56 ±0.02 



dialed 'hot Jupiter' atmospiieres (e.g. HD 189733b. iKnutson et al.l 
l2007b . 

NLTT 5306 (=SDSS J 013532.98+ 1 44555.8) wa s first iden- 
tified as a candida te WD+BD binary in Isteele et all l l201lh and 
iGirven et al l|). The former used an estimate of the WD's at- 
mospheric parameters {T^a — 8083 ± 22 and loggr — 8.08 ± 
0.04; Eisenstein et al. 2006.) in combination with cooling models 
for hydrogen atmosphere (DA) WDs ([Ho lberg & Bergeron 2006; 
Kowalski & Saumod [20o3 ; | Tremblav et al., ,201 1, ; .Bergeron et al.i 
201lh to predict the star's near-infrared photometry. A compar- 
ison was then made with the UKIDSS observations identifying 
NLTT 5306 as having a near-infrared excess consistent with a red 
companion. The SDSS and UKIDSS magnitudes are given in Ta- 
ble [T] Further fitting of the photometry yielded an estimated spec- 
tral type of dL5 for the secondary, with a mass of 58 ± 2Mjup 
at a distance of 60 ± 10 pc. It should be noted that this spectro- 
scopic mass estimate is model dependent, c alculated by interpo- 
lating the Lyon gr oup atmospheric models JChabrier et al.llioool ; 
iBaraffe et alj2002l) given an estimated age for the WD and temper- 
ature for the BD. The system was unresolved with an upper limit 
on the projected separation of < 57 AU between components. 

The structure of this paper is as follows; In Section 2 we de- 
scribe the observations and their reduction. In Section 3 we de- 
scribe the analysis of the data, starting with the optical light curve, 
followed by the optical and NIR spectroscopy of the WD and BD, 
and finally the radial velocity. In Section 4 we discuss the implica- 
tions of the results and state our conclusions. 

2 OBSERVATIONS AND DATA REDUCTION 
2.1 INT Optical Photometry 

NLTT 5306 was observed photometrically for two hours (80 x 90s 
exposures) in the Sloan i'-band on the night of 2009 October 23 
with the Wide Field Camera (WFC) on the 2.5 Isaac Newton Tele- 
scope (INT) in La Palma, Spain. The da ta were reduced using the 
INT Wide Field Survey (WFS) pipeline dlrwin & LewislboOll) de- 
veloped by the Cambridge Astronomical Sur vey Unit. For a de - 
tailed description of the reduction process, see llrwin et sl\ ilOOH) . 
In brief a standard CCD reduction was performed by correcting for 
the bias, trimming the frames, correcting non-linearity, flat-fielding 
and correcting for the gain. The flux was measured in each obser- 
vation using aperture photometry and the result converted to mag- 
nitudes using nigh tly zero-point est imates based on standard star 
field observations dlrwin et al.ll2007h . 



2.2 X-Shooter Spectroscopy 

NLTT 5306 was observed using X-Shooter dD'Odorico etal.l2006h 
mounted at the VLT-UT2 telescope on the night of 2010 Septem- 
ber 5. X-Shooter is a medium-resolution spectrograph capable of 
observing using 3 independent arms simultaneously; the ultraviolet 
(UVB), optical (VIS) and the near-infrared (NIR) arms covering a 
wavelength range of 0.3-2. 5/im. For our observations we used slit 
widths of 0.8", 0.9" and 0.9" in the UVB, VIS and NIR arms re- 
spectively. Exposure times for each arm were 1200 s in the UVB, 
1200 s in the VIS and 12 x 150 s in the NIR. We nodded between 
each exposure along the NIR slit to improve sky subtraction. This 
gave us a total of 4 exposures in each arm. 

Reduction of the raw frames was carried out using the stan- 
dard pipeline release of the ESQ X-Shooter Common Pipeline Li- 
brary recipes (version 1.3.7) within GASGANcQ, version 2.4.0. 
The standard recipes were used with the default settings to reduce 
and wavelength calibrate the 2-dimensional spectrum for each arm. 
The extraction of the science and spectrophotometric standard and 
telluric spectra were carried out using APALL within IRAF. The 
instrumental response was determined by dividing the associated 
standard star by its corresponding flux table. We also used this 
method to apply the telluric correction. 

Finally, the spectra were flux calibrated using the SDSS and 
UKIDSS magnitudes (Table[lll. 

2.3 SDSS Spectroscopy 

NLTT 5306 was observed on multiple occasions by the SDSS 
(=SDSS J013532.97-H144555.9). From the SDSS archive we ex- 
tracted a total of 17 spectra taken over the period from 2000 De- 
cember 1 to 20. All but three of these spectra had exposure times 
of 15 minutes; the three spectra acquired on December 4 were ex- 
posed for 20 minutes. 

2.4 High Resolution Spectroscopy with the Hobby-Eberly 
Telescope 

NLTT 5306 was observed using the High Resolution Spectrograph 
(HRS; Tull 1998) on the Hobby-Eberly Telescope (HET; Ramsey 
et al. 1998) on the nights of 2010 December 6 & 10, 201 1 January 
16 & 23, and 2011 February 8. The ephemeris of the system was 
unknown at the time and so 6 random observations were taken in 
order to establish if the primary has a measurable radial velocity 
variation, and to then estimate an orbital period. Each observation 
was split into 2 separate exposures of 1320 s. A ThAr lamp ex- 
posure was taken both before and after the science observations in 
order to aid wavelength calibration. 

The cross disperser setting was '316g5936' corresponding to 
a wavelength range of 4076-7838A, in order to cover the Ho, H/3 
and Balmer lines. This gives a spectral resolving power of 
R = X/SX — 15, 000. Two sky fibers were used to simultaneously 
record the sky background. 

Reduction of the raw frames was carried out using standard 
routines in IRAF. In brief, bias and flat frames were combined 
and used to correct the science frames. The extraction of the sci- 
ence spectra were carries out using APALL within IRAF. The sky 
spectrum was extracted in the same way as the science, simply by 
shifting all the apertures by a set amount to cover the parallel sky 

^ http://www.eso.org/sci/software/gasgano 
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Table 2. Measured atmospheric parameters for NLTT5306 A from various spectroscopic observations. 



Telescope/Instrument 


Teff (K) 


logs 


M/rrriWD ( Mq) 


d(pc) 


SDSS 1 


7641 ± 48 


7.61 ±0.10 


0.39 ±0.05 


72 ±4 


SDSS2 


7729 ± 7 


7.67 ±0.05 


0.42 ±0.02 


68 ± 2 


SDSS 3 


7729 ± 49 


7.71 ±0.10 


0.44 ± 0.05 


68 ±4 


VLT+X-Siiooter 


7925 ± 15 


7.74 ± 0.02 


0.47 ±0.01 


70 ± 1 
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Figure 1. Floating-mean periodogram for the INT i'-band photome- 
try of NLTT5306. The global minimum is located at a frequency of 
14.115-4 cycles/day. 
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Figure 2. Phase folded INT i'-band light curve of NLTT5306 showing a 
peak-to-peak variability of ~ 1% with a period of 102. 2^ min. 



fiber. The extracted sky spectrum was then scaled so that the sky 
lines matched the strength of the corresponding lines in the science 
spectra. Particular attention was payed to the order containing the 
Hq absorption. This was then subtracted from the science spectra. 
Finally the sky subtracted science spectra were normalised using 
the CONTINUUM package within IRAF. 



3 ANALYSIS 

3.1 Optical Light Curve 

The INT i'-band light curve of NLTT5306 (Figure |2) shows low- 
level photometric variability with a peak-to-peak amplitude of ~ 
0.8%. The analysis leading to this result is outlined below. 

A 'floating- mean' periodogram was used to sea rch for period- 
icity in the target jCumming. Marcv & Butle3ll99 9'). This method 
involves fitting the time-series data with a sinusoid plus a constant 
A in the form of: 

A + Bsin[27r/(t - to)] (1) 

where / is the frequency and t is the time of observation. The re- 
sulting periodogram is a plot of the fit with frequency (Figure[T). 
The errors associated with the best-fitting frequency were estimated 
as the 2a frequency range from the global minimum (correspond- 
ing to a change in of four, assuming only one useful fitted pa- 
rameter). 

To evaluate the significance of the best-fitting period, a false 
alarm probability (FAP) was estimated using 100000 Monte Carlo 
trials and an analytical approach for comparison. Fake light curve 
datasets were generated for the Monte Carlo tests at the same 
timings as the observations with the mean magnitude as the ob- 
served data. Random Gaussian noise was then added to the flux 
distributed with the same variance as the observed magnitudes. The 
FAP was determined from the number of trials where the maxi- 
mum power in the periodogram (from the fake dataset) exceeded 
the maximum power from the observed dataset. A significant de- 
tection threshold was set at 1% (a of FAP ^ 0.01). The analytical 
probability was determ i ned using the equation given in Table 1 in 
IZechmeister & KUrsterl ( l2009j) for the residual vari ance normalisa- 
tion (also see Appendix B in lCumming et al .1 19991). Further details 
on the significance tests can be found in ICumming et"^ ( 1 19991) and 
Lawrie et al. (2012, in prep.). 

A global minimum in the periodogram is found at a frequency 
of 14.lj;?;4 cycles/day, a period of 102.21 "4^ min. A fitted sine 
wave to the data gives a reduced of 2.07 (x^ of 157 over 76 
degrees of freedom (dof)), while a constant fit to the data gives a 
reduced x^ of 2.92 (x^ of 231 over 79 dof). The FAP statistics 
are well within the limit for a significant detection, with a FAP of 
<0.001 from the Monte Cario tests and a FAP of 5 x 10"^ from the 
analytical estimation. This suggests that it is significantly unlikely 
that the variability seen in the light curve is due to noise fluctuations 
alone. 

The probability that this system has an alignment that would 
result in an eclipse as viewable from Earth is ~ 20%, with an 
eclipse duration of ~ 7mins (Faedi et al. 201 1). No eclipse or graz- 
ing transit is immediately obvious in the phase folded light curve 
of NLTT5306 (Figure |2j. Given the total time between exposures 
of 90s and the coverage of ~ 1.5 orbital periods, it seems unlikely 
that we would have missed an eclipse. 
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0.6 0.7 
Wavelength (micron) 

Figure 3. X-Shooter spectram (grey) of WD0132+142 covering the UVB 
and VIS arms where the WD primary dominates. SDSS ugriz photometry 
is overplotted (squares). The WD model SED (black) begins to diverge from 
the observed spectram at longer wavelenghts due to the added flux from the 
secondary. 
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Figure 5. Balmer line fit (black) of NLTT 5306 A using the X-Shooter spec- 
tram (grey). The emission line seen at the core of Ho is likely due to accre- 
tion via a stellar wind. The region of the emission line was excluded from 
the fit. 



0.6 
0.55 

0.5 
0.45 



0.3 
0.25 - 
0.2 




1.4 1.6 1.8 2 
Wavelength (micron) 

Figure 4. X-Shooter Spectrum (grey) of WD 0132+142 covering the NIR 
arm where the BD secondary dominates. UKIDSS YJHK photometry is 
overplotted (squares) as well as a composite WD+dL5 model (black dashed) 
for visual compaiison. 
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Figure 6. X-Shooter spectram of NLTT 5306 B (grey) created by subtract- 
ing the WD model spectram calculated using the measured atmospheric 
parameters of the primary. 



3.2 White Dwarf Spectroscopy 

Figures |3] and |4] show the extracted X-Shooter spectmm for 
NLTT 5306 separated into the UVBA^IS and NIR arms. This 
clearly shows the WD primary dominating the optical wavelengths 
and the BD companion dominating the near-infrared wavelengths. 
A closer inspection of the Ha absorption line in individual expo- 
sures revealed the presence of line emission close to the line centre 
(Figure|5}. 

Figure |5] and Table |2] shows the results of fitting the Balmer 
series with atm ospheric models of hydrogen atmosphere DA WDs 
( lKoesteJl2008l) . for both the SDSS and X-Shooter spectra. The av- 
erage of these values yields and effective temperature and surface 
gravity of T^a = 7756 ± 35K and logg = 7.68 ± 0.08 respec- 
tively. We interpolated these values over agrid of synthetic colours 
and evolutionary sequences of DA WD^ to calculate a mass of 
MwD = 0.44 ± 0.04 Mq and a distance of 71 ± 4pc. The re- 
sults are summarised in Table |4] The region within Ha containing 



the core emission was excluded from this fit. A model spectrum 
calculated using these values is overplotted in Figure|3] 

It should be noted that the spectroscipic fit may have lead to an 
overestimate of the WD's mass, due to the difficulty in modelling 
the hydrogen line profiles below T^s ^ 12000 - 13000 K. The 
iTremblav. Bergeron & GianninasI (1201 ih analysis of SDSS WDs 
suggest an overestimate at the level of 10%, although this is largely 
based on WDs with a mass of around 0.6 Mq. A lower WD mass 
would decrease the minimum value for the mass of the secondary, 
relaxing the need for the system to have a high inclination for the 
mass to be consistent with the spectral type of dL4-7 and given the 
non-detection of an eclipse or grazing transit (Section r3.lt . 



3.3 Brown Dwarf Spectroscopy 

The model WD spectrum for NLTT 5306 A was subtracted from the 
observed X-Shooter spectrum in order to estimate the spectral en- 
ergy distribution of solely the secondary, NLTT 5306 B (Figure|6j. 

In order to estimate a spectral type for the secondary we have 
calculated the H2O'' and H20'^ indices, defined by: 



http://www.astro.umontreal.ca/~bergeron/CoolingModels/ 



H2O" 



F(1.48)/F(1.60) (2) 
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Table 3. Observation times and measured radial velocities from the B aimer 
absorption lines of NLTT5306. The SDSS data have been shifted to the 
mean systemic velocity of the X-Shooter and HET observations. 
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Figure 7. H2O'' and H20'^ Indices as measured for L-dwarfs from the NIR- 
SPEC Infrared Archive (crosses). Ths short-dashed lines shows the same 
indices measured for NLTT 5306 B. These suggests the spectral type in the 
approximate range dL4-dL7. The H2O'' index (Upper panel) is the better 
indicator of the spectral type. Equation 4 is overplotted in the upper panel 
(long-dashed) to show the linear relationship between spectral type and the 
H2O'' index. 



H2O" = F(1.80)/F(1.70) (3) 

where F is the average flux in a band, w i th ±0 .01^m, centred on the 
specified wavelength. iBurgasser et al.l (I2OO2I) state that the H2O'' 
index in particular shows a linear relation with spectral type given 
by: 

SpT = (12.6 ± 0.9) - (26.7 ± 0.6)H2O'' (4) 

where SpT=0 at dTO and SpT=-4 at dL5. This relation holds for 
spectral types between dM5 (-14) and dT8 (8). 

Figure [7] shows these indices as measured for L-dwarfs from 
the NIRSPEC infrared archivfH as well as the H2O'' = 0.637 and 
H2O'' = 0.627 as measured for NLTT5306 B. Substituting the for- 
mer value into equation 4 gives SpT=— 4.4 ± 1.3, corresponding to 
a spectral type between dL4-dL7. This is consistent w ith the origi- 
nal ph otometry based spectral type estimate of dL5 in lSteele et al.l 
( I2OIII) . 



3.4 Radial Velocity 

Each SDSS spectrum is a combination of separate red and blue 
components. We split these into two datasets, the first covering Ha 
and the second covering the higher Balmer series. We fitted mul- 
tiple Gaussian components to these lines, using two Gaussians per 
line (see Marsh, Dhillon & Duck 1995 for details of this process). 
Once profiles to the mean spectra had been fitted, then for the final 
fit we held all shape parameters fixed and simply allowed the radial 
velocity to be fitted, giving us our final radial velocities. 

The HET spectra also consists of separate blue and red com- 
ponents, the blue end covering HjS and H7, and the red covering 



http://www.astro.ucla.edu/ mclean/BDSSarchive/ 



Telescope/Instrument 



HJD 



RV(kms-i) 



SDSS 2451879.73722 


-39.06±14.4 


2451879.74935 


-57.86±14.6 


A.'-rJ io 1 y . 1 oJy I 


f,AA- 1 7 7 


J i OoZ. / JtZ? 


S 1 lAA-Q R 
J i . A^^y . 


/.H-J 1 ooZ,. / / UUO 




2451882.78609 


-46.46±11.5 


2451884.67157 


48.14±I1.0 


2451884.68378 


-3.96±11.7 


2451884.69593 


-52.26±10.9 


2451884.70808 


-17.46±12.2 


2451884.72014 


10.94±12.2 


2451884.73226 


44.64±13.8 


2451884.74498 


37.34±11.9 


2451884.75713 


-12.76±11.9 


2451898.68069 


28.84±8.6 


2451898.69285 


-11.86±9.1 


2451898.70504 


-46.76±9.0 


VLT+X-Shooter 2455444.86040 


-43.01 ±4.50 


2455444.86933 


-12.61±4.94 


2455444.87819 


16.90±4.72 


2455444.88711 


43.71±4.65 


HET+HRS 2455536.74182 


-26.24±4.13 


2455536.75732 


-59.02±6.34 


2455540.72758 


-29.01 ±4.09 


2455540.74309 


19.97±6.01 


2455577.62159 


27.81±4.50 


2455577.63707 


.9.44±4.58 


2455584.61300 


59.44±5.I1 


2455584.62852 


51.29±8.52 


2455600.57543 


-23.57±5.74 



Hq. Four of the five observations consisted of 2 separate exposures 
of 1320s each, with the final observation only producing one such 
exposure. As these exposures sampled a significant fraction of the 
orbital period of NLTT 5306 B, each observation was fitted sepa- 
rately, giving a total of 9 data points. Four more radial velocities 
were measured by extracting individually the spectra which were 
combined to make the final X-Shooter spectrum. 

The absorption lines in the HET and X-Shooter data were fit- 
ted using a combination of 2 Gaussians using the program FITS B 2 
written by Ralf Napiwotzki. The best fitting shape parameters were 
fixed and then fitted for velocity shifts. 

The final measured radial velocities for all 3 datasets are given 
in Table [3] The SDSS data were shifted so that the mean veloc- 
ity matched that of the higher resolution and therefore superior 
HET and X-Shooter data. The low resolution of the SDSS data 
made measurement of the systemtic velocity with just the SDSS 
data unreliable, with a significant difference measured in the mean 
velocites ('^ 20km/s) even between the blue and red arms. The 
mean velocity measured from HET and X-Shooter is consistent 
with the systemic velocity of -15±36 km/s previously measured 
in Adelman-McCarthy (2008). 

The high resolution of HET HRS and X-Shooter spectra al- 
lowed for the detection of the line emission in the core of Ha. The 
fit to the WD's Ha absorption line was used to subtract the contri- 
bution of the WD. The emission line was then fitted with a single 
Gaussian to measure the radial velocity. This was possible in all 4 
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Figure 8. Periodgram covering the region of interest for tlie WD's ra- 
dial velocities as measured using the Balmer series absorption lines in the 
SDSS, HET and X-Shooter data. Two aliases are favoured at periods of 
101.88±0.02 and 109.96±0.02 mins with a minimal difference in ■ 
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Figure 10. Phase folded radial velocity curve for the favoured period of 
101.87±0.04mins of the H absorption lines of NLTT 5306 A showing the 
SDSS (squares), HET (crosses) and X-Shooter (circles) data. The curve rep- 
resenting the best fitting orbital parameters is over-plotted (dashed). 
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Figure 9. Periodgram covering the region of interest for the BD's radial ve- 
locities as measured using the Ha core emission in the HET and X-Shooter 
data. In this case the periodogram favours the higher frequency alias giving 
a period of 101.87±0.04 mins. 
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Figure 11. Phase folded radial velocity curve for the favoured period of 
101.87±0.04 mins of the Ho emission line showing the HET (crosses) and 
X-Shooter (circles) data. The curve representing the best fitting orbital pa- 
rameters is over-plotted (solid) as well as the best fitting model for the ab- 
sorption lines (dashed). 



X-Shooter exposures and 8 of the 9 HET exposures, with the 8th 
exposure being too noisy to accurately identify the emission. 

Periodicity was searched for in the combined radial velocity 
data (Table ID of the WDs' absoiption lines (SDSS, HET, and X- 
Shooter) and for the line emission (HET and X-Shooter) using the 
method discussed in Section lSTI For the WD, 30 unique data points 
were used as measured from the Balmer series absoiption lines, and 
12 as measured from the line emission in the core of Ha. 

Figure [8] shows the periodogram for the fit of the WD radial 
velocity measurements. The analysis of this data resulted in 2 fa- 
vored aliases with periods of 109. 96±0.02 mins (X^ = 72.6) and 
101. 88±0.02 mins (x^ = 66.5). The x^ value slightly favors the 
shorter period alias which would be consistent with the period mea- 
sured from the INT light curve. However, the longer period alias 
produces a period which is also within errors of this value. There- 
fore, it was not possible to determine the period accurately just util- 
ising the available WD data. 

Figure |9] shows the periodogram for the fit of the Ha core 
emission radial velocity measurements. The 2 most favoured pe- 
riods of 101.87±0.04mins (x^ = 78.5) and 109.64±0.04 mins 



(X^ ~ 93.0) are indeed consistent with those obtained for both 
the light curve and the radial velocity fit of the WD. However, the 
X^ value is much more in favour of the shorter period alias in this 
case. Therefore, we take the value of P = 101. 88±0.02 mins from 
the larger (and therefore more accurate) absoiption line dataset for 
the period of the system. 

An attempt was made to measure the radial velocity of the 
BD component using the four individual X-Shooter observations 
in the NIR arm. Each spectrum was extracted and flux calibrated, 
and the telluric correction applied. The model WD spectrum (Fig- 
ure[3} was then subtracted from each spectrum as in Section 3.3. We 
then cross correlated each spectrum with a dL5 template using the 
IRAF package FXCOR. This was attempted using various wave- 
length ranges in order to negate the water vapour bands. However 
the S/N in the individual exposures was such that the cross corre- 
lations with the template did not produce any sensible results. We 
also suspect that we were not entirely able to remove the WD con- 
tribution satisfactorily, and it is also possible that the BD spectrum 
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undergoes short terai changes in overall shape due to the irradiation 
from the primary. 



4 DISCUSSION 

The s tellar/substell ar nature of an object is dependent on its 
mass jKumaill 19631) . The commonly used limit to distinguish be- 
tween low mass main-sequence stars and BDs is O.O75M0(75Mj„p, 
iBurrows et all 1 19971 : IChabrier & BaraffelbOOol) . below which hy- 
drogen fusion does not occur. 

Figure [To] and [TT] show the phase folded radial velocity mea- 
surements for NLTT5306. The parameters of the spectroscopic or- 
bit are summarised in Table|5] The period measured using the radial 
velocity data is consistent with the period measured using the vari- 
ability of the INT i'-band light curve. Using the calculated value 
for the mass of the primary, Mwd =0.44±0.04 M© (Table©, the 
minimum mass of the secondary is 56±3Mjup, consistent with the 
measured spectral type of dL4-dL7 and confirming NLTT 5306 B is 
a bona fide BD. Given this consistency, we suspect that this system 
has a relatively high inclination, even though no eclipse was de- 
tected. It also co nfirms NUTT5306 is the shortest period detached 
WD+BD binarv. lParsons et al.l ( l2bl2al) detected a detached WD bi- 
nary system with a period of only 94mins (CSS 03170). However, 
in this case the secondary has a mass above the hydrogen burning 
limit and so is classed as a main sequence star. The most simi- 
lar known system is WD 0137— 349 which has an orbital period of 
116mins (Maxted et al. 2006) and a slightly lower mass secondary 
(53±6Mj„p, Burleigh et al. 2006). 

One might assume that the Ha emission seen in our spec- 
troscopy is likely to arise from the irradiation of the BDs atmo- 
sphere by the WD primary, as is the case for WD 0137— 349. If 
this were true, then the radial velocity measurements would be in 
anti-phase with those measured from the Ha absorption, and the 
amplitude would allow us to solve for the masses of both binary 
components. However, Figure [TT] shows the emission is clearly in 
phase with the absorption, and with a similar measured amplitude 
of 48.9±1.8 km/s, its origin must be associated with the WD. The 
most likely cause of such emission is accretion, either via Roche 
Lobe overflow or wind from the substellar companion. 

Burleigh et al. (2006b) observed a similar situation in the 
magnetic WD+BD binary SDSS J 121209.31+013627.7 (Porb « 
90mins). In this case the system is considered to be in a semi- 
detached state with a magnetic cataclysmic variable (polar) in a 
low state of accretion from a BD onto a magnetic WD. Although 
NLTT 5306 shows no evidence of a detectable magnetic field and at 
Porb ~ 102 mins is probably not in semi-detached contact, it may 
be accreting from a weak wind from the brown dwarf. Therefore, 
NLTT 5306 is more akin to the wind accreting system LTT560 
( jTappert et ai]l2007l I20T l|), albeit in this system the secondary is 
a much earlier spectral type of dM5. 

The Ha emission seen in LTT 560 consists of two anti-phased 
components, one originating in the secondary and the other from a 
chromosphere on the WD as a result of accretion via the compan- 
ion's stellar wind. NLTT 5306 only shows the emission component 
associated with the WD, and so we conclude we may be observing 
a similar situation where there is only chromospheric emission and 
no obvious activity from the secondaiy (which is to be expected 
given the estimated spectral type). The origin of this emission line 
component in post-common-env elope binaries (PCEB) is briefly 
discussed in lTappert et al.l ( 1201 ll) , although given the rarity of sys- 
tems where the Ha emission line component is located on the WD, 



it is presently unclear under what conditions chomospheric emis- 
sion occurs. 

If NLTT 5306 is accreting via a stellar wind onto a chromo- 
sphere then this would occur some distance above the WD. The 
systemic velocity of the observed emission is similar in value to that 
of the radial velocity measured using the WDs Balmer absorption 
lines. Some difference would be expected due to the gravitational 
redshift of the WD if the emission was chromospheric in origin. 
Adopting the measured atmospheric parameters from this work, 
this amounts to a redshift of Vgr ~ 17.9 ± 3.5kms~^. Figure 11 
shows the systemic velocity of the emission feature is ^ 5 — lOkm/s 
greater than that measured from Balmer absorption. If anything we 
would expect it to be less (i.e redshifted compared to the WD) but 
given the low S/N of the emission detected in the HET data (from 
which the fit heavily relies), this value should not be trusted. In 
all likelihood the emission forming region is somewhere above the 
WDs photosphere. This could be further constrained with more ac- 
curate measurements of the radial velocity of the absorption and 
emission features, and should be considered as a future project for 
the VLT + X-Shooter. 

The effective temperature of the BD can be estimated to be 
~ 1700 K from its measured spe ctral type of dL4 -dL7. and com- 
parison with observed L-dwarfs dVrba et al.ll2004l) . Using this ef- 
fective temperature and the cooling age of the WD (Table |4j as a 
minimum value for the age of the system, we have estimated the 
radius of NLTT 5306 B to be _Rbd = 0.95 db 0.04 Rjup by inter- 
polating the Lyon g roup atmospheric models dChabrier et al.l2OO0l : 
iBaraffe et alj2002i) . 

The systemic velocity of the WD allows us to discuss the kine- 
matics of the system, in particular the U velocity. This gives a good 
indication of whether the WD is a thin disc, thick disc or halo ob- 
ject, and thus allows for fur ther constraints on the age of the pri- 
mary. Using the equations o f jjohnson & SoderblomI fl987h and the 
values given in Tables [4] and [5] we calculate U ~ 70kms~^ for 
NLTT 5306. This would seem to suggest that t he WD is a mem- 
ber of the thick disk population (See Figure 4 of IPauli et al.ll200^ 
and is likely much older than the minimum cooling age suggests 
(> 5Gyr). At this age the Lyon group models give a mass of closer 
to 70 Mjup for the companion, still well within the accepted BD 
range. 

The asymmetric heating and rotation of the BD pro duces a 
modulation of brightness known as the 'reflection effect' (IWilsoiJ 
119901) . This has been detected at the order of ~ 1% in the INT i'- 
band light curve. Since the BD is tidally locked with the WD, this 
has allowed us to estimate the binary orbital period independently 
of the radial velocity measurements. A more accurate spectral typ- 
ing of the companion would require further measurements of this 
effect at longer wavelengths (i.e. t he near-infrared) whe re this vari- 
ation would be more pronounced ( iBurleigh et al.l2008h . 

Variability has also been observed in the i'-band for 
WD 0137-349, but of the order of ~ 2% (Burleigh, pri- 
vate communication). WD0137-349 A is hotter than NLTT 5306 
(~ 16000 K) so there are more UV photons, and more flux overall 
by approximately an order of magnitude. Therefore, we would not 
necessarily expect to observe such a strong effect on NLTT 5306 B, 
which may also depend on local conditions and chemistry in the 
BD atmosphere, but a variation of ~1% seems consistent with the 
effects of irradiation. 

The progenitor system of NLTT5306A&B consisted of a 
main-sequence star and a BD with an orbital separation suffi- 
ciently small for the progenitor of NLTT 5 306 A to fill its Roche 
lobe as it evolved off the main sequence. As a consequence 
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Table 4. Properties of the white dwarf NLTT 5306 A 



Parameter 


Value 


Toff (K) 


7756±35 


logg (c.g.s units) 


7.68±0.08 


Mass (Mq) 


0.44±0.()4 


Cooling Age (Myr) 


710±50 


Radius (Rq) 


().0156±0.0016 


Distance (pc) 


71±4 



of the ensuing unstable mass transfer the BD was engulfed in 
the envelope of the progenitor of NLTT 5306 A, leading to a 
rapid reduction in the orbital period and the ejection of the en- 
velope. The low mass of NLTT 5306 A suggests that the core- 
growth was truncated by this com mon envelope evolution, and that 
this WD may contain a He-core (IWebbinldll984l : llben & Tutukovl 
1 19861 ; iRebassa-Mansergas et al.ll201lh . Therefore, the evolution of 
NLTT 5306 A likely terminated on the red giant branch (RGB) 
rather the the asymptotic giant branch (AGB). 

Following the emergence from the common envelope, the bi- 
nary continued to evolve towards shorter periods. Given the low 
mass of NLTT 5306 B gravitational wave radiation is likely to be 
the only relevant agent of orbital angular momentum loss. Adopt- 
ing the stellar parameters for the WD and the BD determined above , 
and using the formalism outlined bv lSchreiber & Gansickel ( l2003l) . 
we calculate the orbital period at the end of the common envelope to 
have been Pce — 120 min. The orbital period of NLTT 5306 will 
continue to decrease for another ~ 900 Myr, until the BD will even- 
tually fill its Roche lobe and initiate stable mass transfer onto the 
WD. This transformation into a CV will occur at an orbital period 
of ~ 68 min, near the orbital period minimum of CVs (Gansicke et 
al. 2009). 

The exis tence of WDO 1 37-349 B ("Maxted et alj l200el : 
iBurleigh et all ,2006) and NLTT 5306 B demons trate that BDs 
can s urvive common envelope evolution (see also iNordhaus et al] 
l2010l) . and their short orbital periods and low WD masses are in 
line with the statistics of the much larger sample of post-conmnon 
envelo pe binaries containing low-mass M-dwarfs (Zorot ovic et al.l 
l2011ah . Binary population models predict both the existence of 
CVs bom at very short periods with BD donors (Politano 2004b), 
and CVs containing low-mass He-core WDs (e.g. de Kooliil992i : 
|PoUtand [r996). Yet, among the sample of known CVs, there is no 
compelling evidence for ei ther systems that were bom with a BD 
donor (see the discus s ion inlLittlefair et al .120071 : lUthas et alj201ll : 



|Parsonsetal]|2012bl: 



Breedt et al 



2012h or CVs conta ining low- 



mass WDs jZorotovic et al.ll2011bl ; ISavourv et al.ll201 l]). We con- 
clude that WD0137-349 and NLTT5306 represent nearby bona- 
fide progenitors of CVs with low-mass WDs and brown-dwarf 
donors, and that the lack of such systems among the CV popula- 
tion reflects that the present-day population of pre-CVs is not fully 
representative of the progenitors of the present-day population of 
CVs. 



5 SUMMARY 

We have spectroscopically confirmed the shortest period WD+BD 
binary known to date. Radial velocity variations and i'-band vari- 
ability due to 'day' and 'night' side heating of the secondary give 
us a period of 101.88 ± 0.02 mins, and a minimum mass for the 
companion of 56±3A/j„p. This is consistent with the spectral type 



Table 5. Spectroscopic orbit of NLTT 5306 where the WD radial velocity 
at a time T is given by 71 + Ki sin[27r/(T — To)], and the emission 
line radial velocity 72 + K2 sin[27r/{T ~ Tq)], where / = 1/P is the 
frequency. 



Parameter 


Value 


P (mins) 


101.88±0.02 


To (HJD) 


2453740. 1408±0.0005 


Ki (kms-l) 


48.1±1.3 


K2 (kms-l) 


48.9±1.8 


71 (kms"-*^) 


2.74±1.3 


72 (kms~^) 


ll.lOitl.O 


a(RQ) 


0.566±0.005 



estimated from the spectroscopy of dL4-dL7. The results are sum- 
marised in Table |5] Emission near the core of Ha indicates accre- 
tion either via a stellar wind or Roche Lobe overflow. NLTT 5306 B 
has survived a stage of common envelope evolution, much like its 
longer period counterpart WD 0137— 349. Both systems are likely 
to represent bona-fide progenitors of cataclysmic variables with a 
low mass white dwarf and a brown dwarf donor. 
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